Process and system including an optimization engine with evolutionary surrogate-assisted prescriptions

ABSTRACT

A surrogate-assisted evolutionary optimization method, ESP, discovers decision strategies in real-world applications. Based on historical data, a surrogate is learned and used to evaluate candidate policies with minimal exploration cost. Extended into sequential decision making, ESP is highly sample efficient, has low variance, and low regret, making the policies reliable and safe. As an unexpected result, the surrogate also regularizes decision making, making it sometimes possible to discover good policies even when direct evolution fails.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application claims the benefit of priority to U.S. Provisional Patent Application No. 62/824,909, entitled “OPTIMIZATION ENGINE WITH EVOLUTIONARY SURROGATE-ASSISTED PRESCRIPTIONS” which was filed on Mar. 27, 2019 and is incorporated herein by reference in its entirety.

Additionally, one skilled in the art appreciates the scope of the existing art which is assumed to be part of the present disclosure for purposes of supporting various concepts underlying the embodiments described herein. By way of particular example only, prior publications, including academic papers, patents and published patent applications listing one or more of the inventors herein are considered to be within the skill of the art and constitute supporting documentation for the embodiments discussed herein.

FIELD OF THE TECHNOLOGY

The technology disclosed relates generally to using evolutionary computation to optimize decision making strategies. More particularly, surrogate predictor model is trained, e.g., with historical data, and a prescriptor model, decision maker, is evolved against the surrogate predictor model.

BACKGROUND

Many organizations in business, government, education, and health-care now collect significant data about their operations. Such data is transforming decision making in organizations: It is now possible to use machine learning techniques to build predictive models of behaviors of customers, consumers, students, and competitors, and, in principle, make better decisions, i.e. those that lead to more desirable outcomes. However, while prediction is necessary, it is only part of the process. Predictive models do not specify what the optimal decisions actually are. To find a good decision strategy, different approaches are needed.

The main challenge is that optimal strategies are not known, so standard gradient-based machine learning approaches cannot be used. The domains are only partially observable, and decision variables and outcomes often interact nonlinearly: For instance, allocating marketing resources to multiple channels may have a nonlinear cumulative effect, or nutrition and exercise may interact to leverage or undermine the effect of medication in treating an illness. Such interactions make it difficult to utilize linear programming and other traditional optimization approaches from operations research. A discussion regarding the deficiencies of the prior art process can be found in Creative AI Through Evolutionary Computation by Risto Miikkulainen arXiv:1901.03775v2 (22 Feb. 2020) the contents of which is incorporated herein by reference.

Instead, good decision strategies need to be found using search, i.e., by generating strategies, evaluating them, and generating new, hopefully better strategies based on the outcomes. In many domains such search cannot be done in the domain itself. For instance, testing an ineffective marketing strategy or medical treatment could be prohibitively costly. However, given that historical data about past decisions and their outcomes exist, it is possible to do the search using a predictive model as a surrogate to evaluate them. Once good decision strategies have been found using the surrogate, they are tested in the real world.

Even with the surrogate, the problem of finding effective decision strategies is still challenging. Nonlinear interactions may result in deceptive search landscapes, where progress towards good solutions cannot be made through incremental improvement and thus discovering them requires large, simultaneous changes to multiple variables. Decision strategies often require balancing multiple objectives, such as performance and cost, and in practice, generating a number of different trade-offs between them is needed. Consequently, search methods such as reinforcement learning (RL), where a solution is gradually improved through local exploration, do not lend themselves well to searching solution strategies either. Further, the number of variables can be very large, e.g. thousands or even millions as in some manufacturing and logistics problems, making methods such as Kriging and Bayesian optimization ineffective. Moreover, the solution is not a single point but a strategy, i.e. a function that maps input situations to optimal decisions, exacerbating the scale-up problem further.

A first exemplary RL solution which has been proposed in the existing art is model-based RL which aims to build a transition model, embodying the system's dynamics, that estimates the system's next state in time, given current state and actions. The transition model, which is learned as part of the RL process, allows for effective action selection to take place. While these models allow agents to leverage predictions of the future state of their environment, the model-based RL usually requires a prohibitive amount of data for building a reliable transition model while also training an agent and they can be computationally intractable for all but the lowest dimensional domains. Recent descriptions of exemplary model-based RL solutions can be found in at least the following references: Deisenroth, M. and Rasmussen, C. E. 2011. PILCO: A model-based and data-efficient approach to policy search. In Proceedings of the 28th International Conference on Machine Learning (ICML) (ICML'11). 465-472; Ha, D. and Schmidhuber, J. 2018. Recurrent World Models Facilitate Policy Evolution. In Advances in Neural Information Processing Systems 32 (NIPS'18), Curran Associates Inc., Red Hook, NY, USA, 2455-2467; Wahlstrom, N., Schon, T. B., and Deisenroth, M. P. 2015. From pixels to torques: Policy learning with deep dynamical models. arXiv preprint arXiv:1502.02251 (2015), the contents of which are known to those skilled in the art and incorporated herein by reference.

A second exemplary RL solution which has been discussed in the existing art is model free RL. As a representative model-free off-policy method, Deep Q-Networks (DQN) solves the sample efficiency issue by modeling future rewards using action values, also known as Q values. The Q-network is learned based on a replay buffer that collects training data from real-world interactions. Advanced techniques such as double Q-learning and dueling network architectures makes DQN competitive in challenging problems. Recent descriptions of exemplary model free RL solutions can be found in at least the following references: Mnih, V., Kavukcuoglu, K., Silver, D., Rusu, A. A., Veness, J., Bellemare, M. G., Graves, A., Riedmiller, M., Fidjeland, A. K., Ostrovski, G., and others, Human-level control through deep reinforcement learning. Nature 518, 7540 (2015), 529-533; Hasselt, H. V. 2010. Double Q-learning. In Advances in Neural Information Processing Systems 23, J. D. Lafferty, C. K. I. Williams, J. Shawe-Taylor, R. S. Zemel, and A. Culotta (Eds.). Curran Associates, Inc., 2613-2621; and Wang, Z., Schaul, T., Hessel, M., Van Hasselt, H., Lanctot, M., and De Freitas, Dueling Network Architectures for Deep Reinforcement Learning, In Proceedings of the 33rd International Conference on International Conference on Machine Learning (ICML) (ICML'16), Vol. 48. JMLR.org, 1995-2003, the contents of which are known to those skilled in the art and incorporated herein by reference.

Third, on-policy model-free techniques, policy gradient approaches (sometimes referred to as deep RL) leverage developments in deep neural networks to provide a general RL solution. For example, Asynchronous Advantage Actor-Critic (A3C) in particular builds policy critics through an advantage function, which considers both action and state values. And Proximal Policy Optimization (PPO) further makes actor-critic methods more robust with a clipped surrogate objective. Unfortunately, policy gradient techniques are typically quite sensitive to hyperparameter settings, and require overwhelming numbers of samples, due, among other reasons, to the need to train a policy neural network. Descriptions of existing A3C and PPO approaches may be found in: Mnih, V., Badia, A. P., Mirza, M., Graves, A., Lillicrap, T., Harley, T., Silver, D., and Kavukcuoglu, Asynchronous methods for deep reinforcement learning, In Proceedings of the 33rd International Conference on International Conference on Machine Learning (ICML) (ICML'16). 1928-1937 and Schulman, J., Wolski, F., Dhariwal, P., Radford, A., and Klimov, O., Proximal Policy Optimization Algorithms. CoRR abs/1707.06347 (2017), the contents of which are known to those skilled in the art and incorporated herein by reference.

And there have been descriptions in the existing art of using evolutionary approaches with RL. For example, the evolution of policy gradients has been described in the following recent articles: Houthooft, R., Chen, Y., Isola, P., Stadie, B., Wolski, F., Ho, O. J., and Abbeel, P., Evolved policy gradients, In Advances in Neural Information Processing Systems 31, Curran Associates, Inc., 5400-5409 (2018); Khadka et al, Evolution-Guided Policy Gradient in Reinforcement Learning, 32nd Conference on Neural Information Processing Systems (NeurIPS 2018); and Pourchot, et al., CEM-RL: Combining Evolutionary and Gradient-Based Methods for Policy Search, ICLR 2019.

Finally, evolutionary techniques have been used to directly discover policies. But such direct policy evolution can only leverage samples from the real-world, which is costly and risky since many evaluations have to take place in low-performing areas of the policy space. This issue in evolutionary optimization led to the development of surrogate-assisted evolution techniques and surrogate methods have been applied to a wide range of domains, ranging from turbomachinery blade optimization to protein design as is known to those skilled in the relevant art.

Though today's AI applications handle analytics well, the resulting decision-making is still done manually or by using rigid algorithms designed to optimize one or more outcomes. Insights are provided, but the choice of which action to take is not always obvious. The following real-world example highlights the limits of decision-making with current AI. A company has determined, through the use of advanced modeling and projections, that the company's business is at risk of missing its revenue targets in six months. What should be done? Cut costs? Sell at a discount or renegotiate supplier prices? Retarget sales? Or maybe the right choice is to invest in growth or efficiency. Which combination of these choices is best? Given one of these choices, what is the best way to implement it?

Simple extensions to reporting, such as predictive analytics, produce coarse grain advice. Even if your system determines your costs are high, it's still up to you, the Company, to decide what to do about it. Businesses want to maximize revenue, maximize sales and minimize cost and risk, but these objectives are often not aligned, and focusing on one usually comes at the expense of another. Reducing costs by too much can mean losing valuable talent, and over-emphasis on sales can lead to unsustainably low prices. Neither of those options is viable. The best decisions are made when these factors and objectives are correctly balanced. Further, current AI systems have a hard time effectively dealing with changing conditions. Today, decision-making models often become obsolete more quickly than we expect, with rapid changes in habits, culture, fashion, the economy, technology disruption, and legal and regulatory changes. Models based on historical data also cannot withstand what are often called “black swan” events, i.e., unpredictable events that are beyond what is normally expected of a situation and have potentially severe consequences. Even prior art systems that learn approximately a million permutations are still not agile enough to respond to unexpected states that have not been captured in the training set. They fail miserably if rules are changed.

Accordingly, there remains a need in the art for an approach that can be generalized across many contexts to evolve better policies in a human-like fashion. Such a system would need to effectively draw conclusions by continuously forming and updating mental models of the world, creating hypothetical contexts, and then mentally forecasting the implications of different actions based on historical knowledge of previous outcomes and an understanding of current conditions and how they might change over time. The system would then apply the chosen strategy, observe the results, and then modify and improve the models of the real world.

SUMMARY OF EMBODIMENTS

In a first embodiment, a computer-implemented process for developing an optimized prescriptor model for determining optimal decision policy outcomes includes: building a predictor surrogate model historical training data to predict an outcome; feeding the predictor surrogate model in an evolutionary algorithm framework to train a prescriptor model using evolution over multiple generations, wherein subsequent generations are evolved based on results of prior generations until an optimized prescriptor model is determined.

In a second embodiment, a computer-implemented process for finding a decision policy that optimizes a set of outcomes to solve a known policy problem, includes:

-   -   i. training a predictor model based on historical training data         including context information (C), actions (A) performed in an         included context, and historical C, A, outcome (O) data sets (C,         A, O);     -   ii. evolving prescriptor models with the trained predictor model         as a surrogate over a predetermined number of generations;     -   iii. applying a determined elite evolved prescriptor model at a         predetermined interval to a known policy problem to generate new         decision policy data in the format (C, A, O);     -   iv. collecting the new decision policy data (C, A, O) and adding         the new decision policy data (C, A, O) to the historical         training data; and v. repeating i. through iv. until a         predetermined convergence metric is met.

In a third embodiment computer-implemented process for discovering optimal prescriptor agents in a reinforced learning process in a predetermined domain, comprising:

-   -   i. applying pre-selected prescriptor models in the predetermined         domain;     -   ii. collecting Q value data for each time step for each         pre-selected prescriptor model, wherein a Q value is a reward         vector for the step;     -   iii. training a predictor model based on data collected in ii.;     -   iv. evolving prescriptor models with the trained predictor model         as a surrogate over a predetermined number of generations;     -   v. repeating i. through iv. until a predetermined convergence         metric is met to discover the optimal prescriptor agents in the         predetermined domain.

BRIEF SUMMARY OF FIGURES

The invention will be described with respect to specific embodiments thereof, and reference will be made to the drawings, in which:

FIG. 1 is a generalized Evolutionary Surrogate-Assisted Prescription (“ESP”) model in accordance with an embodiment herein;

FIGS. 2(a)-2(d) are process flows and system schematics of the ESP algorithm and exemplary system components in accordance with an embodiment herein;

FIGS. 3(a)-3(d) provide visualizations of the behavior of ESP in a synthetic function approximation domain in accordance with an embodiment herein;

FIGS. 4(a)-4(b) provide exemplary visualizations of performance of ESP in the function approximation domain in accordance with an embodiment herein;

FIGS. 5(a)-5(b) provide exemplary visualizations of performance of ESP in the cart-pole domain in accordance with an embodiment herein;

FIG. 6 provides an exemplary visualization of a surrogate approximation of the fitness landscape in accordance with an embodiment herein;

FIG. 7 provides an exemplary visualization of performance of ESP in the Flappy-Bird domain in accordance with an embodiment herein; and

FIG. 8 is an alternate process flow of the ESP algorithm including human intervention in accordance with an embodiment herein.

DETAILED DESCRIPTION

To empower decision makers with the right information at the right time warrants design and execution of the right set of interventions based on the data and inputs from various sources. In an era of Artificial Intelligence (AI) and now, Evolutionary AI, it is imperative that businesses utilize these enablers to achieve the stated objective in a short span of time. The present embodiments use evolutionary AI. Evolutionary AI is a general name for a group of problem-solving techniques inspired by the theory of biological evolution, such as genetic inheritance, natural selection, recombination, mutation and survival of the fittest. Techniques in this field are used on problems that have too many variables for traditional algorithms to consider and where the approach to solving a particular problem is not well understood or requires creativity. One skilled in the art is aware of the numerous descriptions in the literature addressing evolutionary AI in its many forms and variations.

There is now significant historical data available on decision making in organizations, consisting of the decision problem, what decisions were made, and how desirable the outcomes were. Using this historical data, it is possible to learn a surrogate model, and with that model, evolve a decision strategy that optimizes the outcomes. The present embodiments introduce a general such approach, called Evolutionary Surrogate-Assisted Prescription, or ESP. The surrogate is, for example, a random forest or a neural network trained with gradient descent, and the strategy is a neural network that is evolved to maximize the predictions of the surrogate model. As discussed further below, ESP can be extended to sequential decision-making tasks, which makes it possible to evaluate the framework in reinforcement learning (RL) benchmarks. Because the majority of evaluations are done on the surrogate, ESP is more sample efficient, has lower variance, and lower regret than standard RL approaches. Surprisingly, ESP solutions are also better because both the surrogate and the strategy network regularize the decision making behavior. ESP thus introduces a foundation to decision optimization in real-world problems.

Through the use of an ESP service within an optimization engine, business outcomes may be impacted by enhancing decision processes. The Learning Evolutionary Algorithm Framework (LEAF) AI-enables manual decision processes using a prescriptive AI system to create and iteratively enhance recommendations, improving processes and achieving business goals in a principled AI-based manner.

At the base of how ESP works are the concepts of Context (C), Actions (A), and Outcomes (O) exemplified in the generalized ESP model shown in FIG. 1 . The process starts with the building of a predictor surrogate model. The predictor can be any machine learning model trained with supervised methods, such as a random forest, a neural network or other models meeting the criteria known to those skilled in the art. The predictor is then used as a surrogate in order to evolve a prescriptor, i.e. a neural network implementing a decision policy that results in the best possible outcomes. The majority of evaluations are done on the surrogate, making the process highly sample-efficient and robust, and leading to decision policies that are regularized and therefore generalize well.

The goal of the ESP approach is to find a decision policy that optimizes a set of outcomes (FIG. 1 ). Given a set of possible contexts (or states) C and possible actions A, a decision policy D returns a set of actions A to be performed in each context C: D(C)=A,  (1) where C∈

and A∈

. For each such (C, A) pair there is a set of outcomes O(C, A), i.e. the consequences of carrying out decision A in context C. For instance, the context might be a description of a patient, actions might be medications, and outcomes might be health, side effects, and costs. In the following, higher values of O are assumed to be better for simplicity.

In ESP, two models are employed: a Predictor P_(d), and a Prescriptor P_(s). As shown in FIG. 2(a), the Predictor takes, as its input, context information (C), as well as actions (A) performed in that context, along with historical C, A, O data sets. The output of the Predictor is the resulting outcomes when the given actions are applied in the given context. The Predictor is therefore defined as P _(d)(C,A)=O′,  (2) such that Σj L(Ij, Oj′) across all dimensions j of O is minimized. The function L can be any of the usual loss functions used in machine learning, such as cross-entropy or mean-squared-error, and the model P_(d) itself can be any supervised machine learning model such as a neural network or a random forest.

As shown in FIG. 2(b), the Prescriptor takes a given context as input, and outputs a set of actions: P _(s)(C)=A,  (3) such that Σi, j Oj′(Ci, Ai) over all possible contexts i is maximized. It thus approximates the optimal decision policy for the problem. Note that the optimal actions A are not known, and must therefore be found through search. The ESP algorithm then operates as an outer loop in a continuous learning and optimization process that constructs the predictor and prescriptor models as follows (FIG. 2 c ):

(1) Train a Predictor based on historical training data; (S1)

(2) Evolve Prescriptors with the Predictor as the surrogate; (S2)

(3) Apply the best Prescriptor in the real world; (S3)

(4) Collect the new data and add to the training set; (S4)

(5) Repeat until convergence. (S5)

As usual in evolutionary search, the process terminates when a satisfactory level of outcomes has been reached, or no more progress can be made. Note that in S1, if no historical decision data exists initially, a random Predictor can be used. Also note that not all data needs to be accumulated for training each iteration. In domains where the underlying relationships between variables might change over time, it might be advisable to selectively ignore samples from the older data as more data is added to the training set in S4. It is thus possible to bias the training set towards more recent experiences.

One skilled in the art will appreciate the myriad of system components, including one or more processing units and databases, and configurations which may be employed to implement the processes and subprocesses described herein. FIG. 2 d is a schematic of an exemplary system. The overall exemplary system includes a surrogate predictor model generation subsystem 10 which includes at base, a database 20 for storing (C, A, O) training datasets, and at least one module for implementing the selected machine learning algorithm 25, which outputs a trained surrogate predictor model 30. The trained surrogate predictor model is used to evaluate fitness of evolved prescriptor model candidates as part of the evolution process implemented by the prescriptor evolution subsystem 40. The prescriptor evolution subsystem 40 includes a prescriptor model candidate generator 45 and a prescriptor candidate population database 50 which can be continually update in accordance with evolving candidates. The prescriptor candidates are evaluated for fitness against the surrogate predictor model 30 by testing module 60 and ranked or otherwise filtered and compared to one another in accordance with the requirements of a competition module 65. Elite prescriptor model(s) 70 are selected for application to real world scenarios by the real world application subsystem 80. A procreation module 55 is used to re-seed and update the prescriptor candidate population database 50 in accordance with known procreation processes. Finally, the outcomes from application of the elite prescriptor model 70 actions to real work scenarios are stored in outcome database 85 and shared with database 20 to update the (C, A, O) training data.

As is appreciated by those skilled in the art, additional modules and processors, servers and databases may be incorporated to perform different or additional tasks, including data processing/filtering/translation as required for different domains and data sources. Further, aspects of the overall system or subsystems may be performed by different entities. For example, the surrogate predictor model generation subsystem 10 and the prescriptor evolution subsystem 40 may be operated by a service provider and provided as a SaaS product, while the real world application subsystem 80 may be operated exclusively by a customer, thus protecting confidential business and other data. The following co-owned patent applications are incorporated herein by reference herein: U.S. patent application Ser. No. 16/424,686 entitled Systems And Methods For Providing Secure Evolution As A Service and U.S. patent application Ser. No. 16/502,439 entitled Systems And Methods For Providing Data-Driven Evolution Of Arbitrary Data Structures.

Building the predictor model is straightforward given a (C, A, O) dataset. The choice of algorithm depends on the domain, i.e., how much data there is, whether it is continuous or discrete, structured or unstructured. In the examples herein, random forests and neural networks will be demonstrated, but one skilled the art will recognize that algorithm choice varies with domain and can select accordingly. The prescriptor model, in contrast, is built using neuroevolution in ESP. Neural networks can express complex nonlinear mappings naturally, and evolution is an efficient way to discover such mappings, and is naturally suited to optimize multiple objectives. Because it is evolved with the predictor, the prescriptor is not restricted by a finite training dataset, or limited opportunities to evaluate in the real world. Instead, the predictor serves as a fitness function, and it can be queried frequently and efficiently. In a multiobjective setting, ESP produces multiple prescriptors, selected from the Pareto front of the multiobjective neuroevolution run.

Applying the ESP framework to RL problems involves extending the contexts and actions to sequences. The prescriptor can be seen as an RL agent, taking the current context as input, and deciding what actions to perform in each time step. The output of the predictor, O′, can be seen as the reward vector for that step, i.e. as Q values (with a given discount factor, such as 0.9, as in the experiments below). Evolution thus aims to maximize the predicted reward, or minimize the regret, throughout the sequence.

The outer loop of ESP changes slightly because in RL there is no dataset to train the predictor; instead, the data needs to be generated by applying the current prescriptors to the domain. An elite set of several good prescriptors are used in this role to create a more diverse training set. The initial training set is created randomly. The loop now is:

(1_(RL)) Apply the elite Prescriptors in the actual domain;

(2_(RL)) Collect Q values for each time step for each Prescriptor;

(3_(RL)) Train a Predictor based on data collected in Step 2_(RL);

(4_(RL)) Evolve Prescriptors with the Predictor as the surrogate;

(5_(RL)) Repeat until convergence.

The evolution of Prescriptors continues in each iteration of this loop from where it left off in previous iteration. In addition, the system keeps track of the best Prescriptor so far, as evaluated in the actual domain, and makes sure it stays in the parent population during evolution. This process discovers good Prescriptor agents efficiently, as will be described in the experiments that follow.

Evolution, which herein refers to evolutionary artificial intelligence technology or evolutionary computation, is a population-based approach. That is, rather than modifying a single solution until it is in an acceptable state (e.g., back-propagating deep networks), the evolutionary computation process in its most general form: generates a population of candidate solutions (initially randomly or using a priori data where available); calculates the fitness of each candidate, removes less fit candidates from the population (e.g., candidate pool) and regenerates new candidate solutions by either randomly tweaking the fitter ones or treating fitter ones as parents and borrowing traits from them and generating new candidates. Thus, evolutionary computation generates a population of models that map a given context (C) to proposed actions (A), and rates the models based on how well they do on the outcomes (O) predicted by the predictor (Pa). Additional features, variations and applications of/to evolutionary computation are described in numerous publications known to those skilled in the art, including at least the following commonly assigned patents and patent applications which are incorporated herein by reference in their entireties: U.S. Pat. No. 8,909,570 entitled Data mining technique with experience-layered gene pool; U.S. Pat. No. 8,918,349 entitled Distributed network for performing complex algorithms; U.S. Pat. No. 8,977,581 entitled Data mining technique with diversity promotion; U.S. Pat. No. 9,002,759 entitled Data mining technique with maintenance of fitness history; U.S. Pat. No. 10,268,953 entitled Data mining technique with maintenance of ancestry counts; U.S. Pat. No. 10,430,709 entitled Data mining technique with distributed novelty search. The above-identified list is not exhaustive.

To test the efficacy of ESP, it was evaluated in three domains: Function approximation, where its behavior could be visualized concretely; Cart-pole control where its performance could be compared to standard RL methods in a standard RL benchmark task; and Flappy Bird, where the regularization effect of the surrogate could be demonstrated most clearly.

The neuroevolution algorithm for discovering Prescriptors evolves weights of neural networks with fixed topologies. Unless otherwise specified, all experiments use the following default setup for evolution: candidates have a single hidden layer with bias and tan h activation; the initial population uses orthogonal initialization of layer weights with a mean of 0 and a standard deviation of 1; the population size is 100; the top 10% of the population is carried over as elites; parents are selected by tournament selection of the top 20% of candidates; recombination is performed by uniform crossover at the weight-level; there is a 10% probability of multi-plying each weight by a mutation factor, where mutation factors are drawn from N (1, 0.1).

The behavior of ESP can be visualized in a synthetic function approximation domain. The domain also allows comparing ESP to direct evolution in the domain, as well as to PPO, and visualizing their differences.

The domain has a one-dimensional context C and a one-dimensional action A, with outcome O given by the function

$O = {- \left| {A - {3\sin\frac{c}{2}}} \middle| . \right.}$ The optimal action for each context lies on a periodic curve, which captures complexity that can arise from periodic variables such as time of day or time of year. The outcome of each action decreases linearly as the action moves away from the optimal action. Episodes in this domain consist of single action in [−10, 10], which is taken in a context drawn uniformly over [−10, 10]. The full domain ground truth is shown in FIG. 3(a).

ESP begins by taking ten random actions. There-after, every iteration, ESP trains a neural network with two hidden layers of size 64 and tan h activation for 2000 epochs using the Adam optimizer with default parameters to minimize MSE, and evolves Prescriptors for 20 generations against this Predictor. Then, the top Prescriptor is run in the real domain for a single episode. Prescriptors have a single hidden layer of size 32 with tan h activation; default parameters are used for evolution.

Direct Evolution (DE) was run as a baseline comparison for ESP. It consists of the exact same evolution process, except that it is run directly against the real function instead of the Predictor. That is, in each generation, all 100 candidates are evaluated on one episode from the real function.

PPO was run as an RL comparison since it is a state-of-the-art RL approach for continuous action spaces. During each iteration it was run for ten episodes since this setting was found to perform best during hyperparameter search. PPO defaults were used for the remaining hyperparameters. Descriptions of PPO which is known to those skilled in the art may be found in Schulman et al., Proximal Policy Optimization Algorithms, arXiv:1707.06347v2 (August, 2017), which is incorporated herein by reference in its entirety.

Ten independent runs were performed for each algorithm. The returned policy at any time for DE and ESP is the candidate with the highest fitness in the population; for PPO it is the learned policy run without stochastic exploration.

Snapshots of the convergence behavior for each method are shown in FIGS. 3(b)-3(d). After 1000 episodes, neither DE nor PPO converged near the optimal solution (FIGS. 3(b)-3(c)). On the other hand, per FIG. 3(d), ESP discovered the periodic nature of the problem within 50 episodes, and converged almost exactly to the optimal within 125 episodes. FIG. 3(d) illustrates ESP Prescriptor (orange) and Predictor (background) for several iterations. The translucent circles indicate the state-action pairs sampled so far, i.e., the samples on which the Predictor is trained. By 125 episodes, ESP has converged around the optimal Prescriptor, and the ESP Predictor has converged in the neighbor-hood of this optimum, showing how ESP can leverage Predictors over time to find good actions quickly. Note that the Prescriptor does not exactly match the actions the Predictor would suggest as optimal: the Prescriptor regularizes the Predictor's overfitting by implicitly ensembling the Predictors evolved over time.

The Predictor's predicted reward for each state-action pair is shown using the background colors in each snapshot of ESP. The rapid convergence of the Predictor highlights the sample efficiency of ESP, due to aggressive use of historical data (shown as translucent circles). Note, however, that the Predictor does not converge to the ground truth over the entire domain; it does so just in the neighborhood of the optimal Prescriptor. Thus, ESP avoids excessive costly exploration of low-quality actions once the structure of optimal actions has become clear.

Note also that the Prescriptor does not follow the optimal action suggested by the Predictor at every iteration exactly. Since it maps states directly to actions, the Prescriptor provides a smoothing regularization in action space that can overcome Predictor overfitting. Also, since the top ESP Prescriptors must survive across many different Predictors over time, ESP benefits from an implicit temporal ensembling, which further improves regularization.

The numerical performance results in FIGS. 4(a) and 4(b) confirm the substantial advantage of ESP. The horizontal axis indicates the total number of real-world episodes used by training and the vertical axis the performance at that point. Ten independent runs were performed for each method. Solid lines represent the mean over 10 runs, and colored areas show the corresponding standard deviation. The true performance of the returned best agent converges substantially faster with ESP (FIG. 4(a)). And ESP also operates with much lower regret than DE or PPO, converging to very low regret behavior orders-of-magnitude faster than the other approaches (FIG. 4(b)). On both metrics, ESP converges orders-of-magnitude faster than the other approaches. In particular, after a few hundred episodes, ESP reaches a solution that is significantly better than any found by DE or PPO, even after 3,000 episodes. This result shows that, beyond being more sample efficient, by systematically exploiting historical data, ESP is able to find solutions that direct evolution or policy gradient search cannot. The standard deviation of ESP is small in both metrics, attesting to the reliability of the method.

Cart-pole experiments described below demonstrate ESP's performance compared to direct evolution and standard RL methods. The Cart-pole control domain is one of the standard RL benchmarks. In the popular CartPole-v0 implementation on the OpenAI Gym platform used in the experiments, there is a single pole on a cart that moves left and right depending on the force applied to it. A reward is given for each time step that the pole stays near vertical and the cart stays near the center of the track; otherwise the episode ends.

DE was run with a population size of 50 candidates. A candidate is a neural network with four inputs (observations), one hidden layer of 32 units with tan h activation, and two outputs (actions) with argmax activation functions. The fitness of each candidate is the average reward over five episodes in the game, where the maximum episode length is 200 time steps.

ESP runs similarly to DE, except that the fitness of each candidate is evaluated against the Predictor instead of the game. A Predictor is a standard multilayer perceptron neural network with six inputs (four observations and two actions), two hidden layers with 64 units each and tan h activation, and one output (the predicted discounted future reward) with tan h activation. It is trained for 1,000 epochs with the Adam optimizer with MSE loss and batch size of 256.

The first Predictor is trained on samples collected from five random agents playing five episodes each. Random agents choose a uniform random action at each time step. A sample corresponds to a time step in the game and comprises four observations, two actions, and the discounted future reward. Reward is +1 on each time step, except for the last one where it is adjusted to +2,000 in case of success, −2,000 in case of failure (i.e. 10× max time steps). The discount factor is set to 0.9. The reward value is then scaled to lie between −1 and 1.

In order to be evaluated against the Predictor, a Prescriptor candidate has to prescribe an action for each observation vector from the collected samples. The action is then concatenated with the observation vector and passed to the Predictor to get the predicted future reward. The fitness of the candidate is the average of the predicted future rewards.

Every five generations, data is collected from the game from the five elites, for five episodes each. The new data is aggregated into the training set and a new Predictor is trained. The generation's candidates are then evaluated on the new Predictor with the new training data. The top elite candidate is also evaluated for 100 episodes on the game for reporting purposes only. Evolution is stopped after 160 generations, which corresponds to 800 episodes played from the game, or once an elite receives an average reward of 200 on five episodes.

In addition to DE and ESP, two state-of-the-art RL methods were implemented for comparison: double DQN with dueling network architectures and actor-critic style PPO. The implementation and parametric setup of DQN and PPO were based on OpenAI Baselines. For PPO, the policy's update frequency was set to 20, which was found to be optimal during hyperparameter search. All other parametric setups of DQN and PPO utilized default setups as recommended in OpenAI Baselines.

FIG. 5(a) shows how the true performance of the best policy returned by ESP, DE, PPO, and DQN changes during the learning process in CartPole-v0. For ESP and DE, the elite candidate that has the best real-world fitness is selected as the best policy so far. For DQN and PPO, whenever the moving average reward of the past 100 episodes of training is increased, the best policy will be updated using the most recent policy. One hundred additional real-world episodes were used to evaluate the best policies (these evaluations are not part of the training).

ESP converges significantly faster than the other methods, implying better sample-efficiency during learning. Moreover, the variance of the true performance becomes significantly smaller for ESP after an early stage, while all other algorithms have high variances even during later stages of learning. This observation demonstrates that the solutions delivered by ESP are highly reliable.

FIG. 5(b) shows the average regret for training processes of all algorithms in CartPole-v0. ESP has significantly lower regret during the entire learning process, indicating not only lower costs but also better safety in real-world interactions.

Whenever a surrogate is used to approximate a fitness function, there is a risk that the surrogate introduces false optima and misleads the search. ESP mitigates that risk by alternating between actual domain evaluations and the surrogate. However, the opposite effect is also possible: FIG. 6 shows how the surrogate may form a more regularized version of the fitness than the real world, and thereby make it easier to learn policies that generalize well.

Flappy Bird is a side-scroller game where the player controls a bird, attempting to fly it between columns of pipes without hitting them by performing flapping actions at carefully chosen times. This experiment is based on a PyGame implementation of this game, running at a speed of 30 frames per second. The goal of the game is to finish ten episodes of two minutes, or 3,600 frames each, through random courses of pipes. A reward is given for each frame where the bird does not collide with the boundaries or the pipes; otherwise the episode ends. The score of each candidate is the average reward over the ten episodes.

Both DE and ESP were setup in a similar way as in the preceding sections. DE had a population of 100 candidates, each a neural network with eight inputs (observations), one hidden layer of 128 nodes with tan h activation, and two outputs (actions) with argmax activation. The ESP Predictor was a random forest with 100 estimators, approximating reward values for state-action pairs frame by frame. The state-action pairs were collected with the ten best candidates of each generation running ten episodes on the actual game, for the total of a hundred episodes per generation.

FIG. 7 shows how the true performance of the best policy returned by ESP and DE improved during the learning process in the flappy bird domain. The elite candidate that has the best real-world fitness was selected as the best policy so far. In about 80,000 episodes, ESP discovered a policy that solved the task, i.e. was able to guide the bird through the entire course of pipes without hitting any of them. It is interesting that DE converged to a suboptimal policy even though it was run an order of magnitude longer. This result is likely due to the regularization effect illustrated in FIG. 6 . Direct evolution overfits to the nonlinear effects in the game, whereas the surrogate helps smooth the search landscape, thereby leading evolution to policies that perform better.

The results disclosed herein show that the ESP approach performs well in sequential decision making tasks like those commonly used as benchmarks for RL. Compared to direct evolution and state-of-the-art RL methods, ESP is highly sample efficient, which is especially important in domains where exploring with the real world is costly. ESP solutions are also reliable and safe, and the complexity of its models can be adjusted depending on the complexity of the data and task. In ESP, only elite agents selected via the surrogate model are evaluated on the real world, significantly improving safety. Second, the quality of the best-recognized policy is unreliable because it has not been sufficiently evaluated during learning. ESP solves this issue by evaluating all elite policies in the real world for multiple episodes. Third, existing RL methods rely heavily on deep neural networks. In contrast, ESP treats the Predictor as a black box, allowing high flexibility in model choices, including simpler models such as random forests that are sufficient in many cases.

These advantages apply to ESP in general, including decision strategies that are not sequential, which suggests that ESP is a good candidate for improving decision making in real-world applications, including those in business, government, education, and healthcare.

When ESP is applied to such practical problems, the process outlined above can be extended further in several ways. First, ESP can be most naturally deployed to augment human decision making. The Presciptor's output is thus taken as advice, and the human decision maker can modify the actions before applying them. FIG. 8 illustrates an exemplary process flow wherein human actions are part of the overall process. These actions and their eventual outcomes are still captured and processed in S4 of the ESP process, and thus become part of the learning (FIG. 2 c ). Second, to support human decision making, an uncertainty estimation model such as RIO can be applied to the Predictor, providing confidence intervals around the out-come O′. Third, the continual new data collection in the outer loop makes it possible to extend ESP to uncertain environments and to dynamic optimization, where the objective function changes over time. By giving higher priority to new examples, the Predictor can be trained to track such changing objectives. Fourth, in some domains, such as those in financial services and healthcare industries that are strongly regulated, it may be necessary to justify the actions explicitly. Rather than evolving a Prescriptor as a neural network, it may be possible to evolve rule-set representations for this role, thus making the decision policy explainable. Such ex-tensions build upon the versatility of the ESP framework, and make it possible to incorporate the demands of real-world applications.

Although the Predictor does not have to be perfect, and its approximate performance can even lead to regularization as was discussed above, it is sometimes the bottleneck in building an application of ESP. In the non-sequential case, the training data may not be sufficiently complete, and in the sequential case, it may be difficult to create episodes that run to successful conclusion early in the training. While the embodiments herein train the Predictor with targets from discounted rewards over time, an alternative approach would be to incrementally extend the time horizon of Predictors by training them iteratively. Such an approach could help resolve conflicts between Q targets, and thereby help in early training. Another approach would be to make the rewards more incremental, or evolve them using reward function search. An yet another approach is to evaluate the quality of the Predictor directly, and adjust sampling from the real world accordingly.

ESP is applicable in multiobjective domains. In many real-world decision-making domains, there are at least two conflicting objectives: performance and cost. As an evolutionary approach, ESP lends itself well to optimizing multiple objectives. The population forms a Pareto front, and multiple Prescriptors can be evolved to represent the different tradeoffs.

ESP is a surrogate-assisted evolutionary optimization method designed specifically for discovering decision strategies in real-world applications. Based on historical data, a surrogate is learned and used to evaluate candidate policies with minimal exploration cost. Extended into sequential decision making, ESP is highly sample efficient, has low variance, and low regret, making the policies reliable and safe. As an unexpected result, the surrogate also regularizes decision making, making it sometimes possible to discover good policies even when direct evolution fails. ESP is therefore a promising approach to improving decision making in many real world applications where historical data is available.

By way of non-limiting examples, the following real-world problems may be addressed using ESP: price and incentive optimization; underwriting; web-site optimization; market-mix optimization; subrogation; fleet management; fraud prevention; human resource management, e.g., hiring, on-boarding, performance, training/development, transitions, wherein outcome objectives such as quality, cost, productivity, time-to-hire, time to performance, post on-boarding overhead, performance and retention are weighed; collections, agricultural growth recipe optimization, e.g., computer-controlled hydroponic growth (objectives: light, temperature, water, nutrients), build surrogate model of how plants grow given a recipe, evolve recipe where each new test is used to improve the model, discover size/taste tradeoff, 24-hour light period; optimizing CRISPR guide RNA; refinery optimization; defect detection and the like.

By way of particular example, the underwriting (C, A, O) data sets might include for Contexts, e.g., business class (e.g., office, store), risk category (e.g., high, low), loss history (e.g., by coverage, location, agent, policy type), policy type (new or renewal), total insured value percentage %; for Actions, e.g., underwriter decision (yes, no); coverage (as-is, new); premium range; and for Outcomes, e.g., insured decision (win, loss), loss ratio (quarterly projection), cost.

Further, elements of the ESP embodiments described herein have been applied for optimizing designs of web pages, where a neural network was used to map user descriptions to most effective web-page designs. Exemplary descriptions are found in Miikkulainen et al., Ascend by Evolv: AI-Based Massively Multivariate Conversion Rate Optimization, AI Magazine (2019) and Miikkulainen et al., Sentient ascend: Ai-based massively multivariate conversion rate optimization, In Proceedings of the Thirtieth Innovative Applications of Artificial Intelligence Conference. AAAI (2018), the contents of which are incorporated herein by reference in their entireties. Similarly, elements of the ESP embodiments described herein have been applied to find effective growth recipes for basil through search with a surrogate model trained with outcomes of past recipes. A description is found Johnson et al., “Flavor-Cyber-Agriculture: Optimization of plant metabolites in an open-source control environment through surrogate modeling,” PLOS ONE (2019), https://doi.org/10.1371/journal.pone.0213918, the contents of which is incorporated herein by reference.

The embodiments described herein combine elements of strategy search and surrogate modeling into a general ESP approach for decision strategy optimization. The embodiments also extend ESP into decision strategies that consist of sequences of decisions. This extension makes it possible to evaluate ESP against other methods in RL domains. Conversely, ESP is used to formalize RL as surrogate-assisted, population based search. This approach is particularly compelling in domains where real-world evaluations are costly. ESP improves upon traditional RL in several ways: It converges faster given the same number of episodes, indicating better sample-efficiency; it has lower variance for best policy performance, indicating better reliability of delivered solutions; and it has lower regret, indicating lower costs and better safety during training. Surprisingly, optimizing against the surrogate also has a regularization effect: the solutions are sometimes more general and thus perform better than solutions discovered in the domain itself. Further, ESP brings the advantages of population-based search out-lined above to RL, i.e. enhanced exploration, multiobjectivity, and scale-up to high-dimensional search spaces.

ESP is a suitable approach to discovering decision strategies when (1) sufficient and consistent historical data of quality measurements on different actions performed in various context is available to bootstrap the surrogate model, (2) deployment allows for reasonable iterations of deploying latest prescriptors and resulting data collection for improving surrogate models, where decision loops currently operated by humans or rigid preprogrammed systems can be replaced.

This application cross-references the following patent applications and publications which are incorporated herein by reference in their entireties and are intended to be part of the present disclosure: U.S. Patent Publication No. 2017/0293849 entitled DISTRIBUTED RULE-BASED PROBABILISTIC TIME-SERIES CLASSIFIER; U.S. Provisional Patent Application No. 62/851,782 entitled Quantifying the Predictive Uncertainty of Neural Networks Via Residual Estimation With I/O Kernel; Qiu, X. et al., Quantifying Point-Prediction Uncertainty in Neural Networks via Residual Estimation with an I/O Kernel, In Proceedings of the Eighth International Conference on Learning Representations (ICLR) (2020); Hodjat, B. et al., PRETSL: Distributed Probabilistic Rule Evolution for Time-Series Classification, In Genetic Programming Theory and Practice XIV. Springer, 139-148 (2018); Meyerson, E. et la., Discovering evolutionary stepping stones through behavior domination, In Proceedings of the Genetic and Evolutionary Computation Conference (GECCO 2017); Miikkulainen, R. et al., Sentient ascend: Ai-based massively multivariate conversion rate optimization. In Proceedings of the Thirtieth Innovative Applications of Artificial Intelligence Conference. AAAI (2018); Miikkulainen et al., Ascend by Evolv: AI-Based Massively Multivariate Conversion Rate Optimization, AI Magazine (2019); Johnson et al., “Flavor-Cyber-Agriculture: Optimization of plant metabolites in an open-source control environment through surrogate modeling,” PLOS ONE (2019), https://doi.org/10.1371/journal.pone.0213918; Stanley, K. et al, “Designing neural networks through neuroevolution” Nature Machine Intelligence, Vol. 1, p.24-35 (January 2019). 

The invention claimed is:
 1. A process for developing an optimized prescriptor model for determining optimal decision policy outcomes comprising: building a predictor surrogate model with historical training data to predict an outcome; and wherein historical training data trains a predictor model, wherein predictor model P_(d) is defined Pd (C, A)=O′, such that Σj L(Oj , Oj ‘) across all dimensions j of O is minimized and function L is a loss function, wherein the loss function (L) is selected from the group consisting of cross-entropy and mean-squared-error, feeding the predictor surrogate model in an evolutionary algorithm framework to train a prescriptor model using evolution over multiple generations, wherein subsequent generations are evolved based on results of prior generations until an optimized prescriptor model is determined, wherein the optimized prescriptor model is determined based on a decision policy, wherein the prescriptor model Ps takes a given context as input, and outputs a set of actions: P_(s) (C)=A, such that Σi, j Oj′ (Ci, Ai) over all possible contexts i is maximized, wherein evolving the prescriptor models over the predetermined number of generations further comprises discovering optimal prescriptor agents in a reinforced learning process in a predetermined domain by; i. applying pre-selected prescriptor models in the predetermined domain; ii. collecting Q value data for each time step for each pre-selected prescriptor model, wherein a Q value is a reward vector for the step; iii. training a predictor model based on data collected in ii.; iv. evolving prescriptor models with the trained predictor model as a surrogate over a predetermined number of generations; v. repeating i. through iv. until a predetermined convergence metric is met to discover the optimal prescriptor agents in the predetermined domain, wherein the pre-selected prescriptor models are neural networks.
 2. The process of claim 1, wherein the historical training data includes context information (C), actions (A) performed in an included context, and historical C, A, outcome (O) data sets (C, A, O).
 3. The process of claim 1, wherein the predictor surrogate model is a machine learning model trained with supervised methods.
 4. The process of claim 3, wherein the predictor surrogate model is a neural network.
 5. The process of claim 1, wherein an elite prescriptor model is iteratively selected and applied to a known policy problem to generate new decision policy data in the format (C, A, O) and the new decision policy data (C, A, O) is iteratively supplied as input to the predictor surrogate model.
 6. The process of claim 5, wherein the elite prescriptor model is selected based on a fitness evaluation.
 7. The process of claim 1, wherein the determination process of the decision policy that optimizes a set of outcomes to solve a known policy problem, comprises: i. training a predictor model based on historical training data including context information (C), actions (A) performed in an included context, and historical C, A, outcome (O) data sets (C, A, O); ii. evolving prescriptor models with the trained predictor model as a surrogate over a predetermined number of generations; iii. applying a determined elite evolved prescriptor model at a predetermined interval to a known policy problem to generate new decision policy data in the format (C, A, O); iv. collecting the new decision policy data (C, A, O) and adding the new decision policy data (C, A, O) to the historical training data; v. repeating i. through iv. until a predetermined convergence metric is met.
 8. The process of claim 7, wherein the predictor model is a machine learning model trained with supervised methods.
 9. The process of claim 8, wherein the predictor model is selected from the group consisting of a neural network and a random forest.
 10. The process of claim 7, wherein the elite evolved prescriptor model is selected based on a fitness evaluation. 